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BRIEF  DESCRIPTION  OF  PORTFOLIO: 

Multi-disciplinary  (Propulsion,  Materials ,  Plasma  and  Electro-Energetic 
Physics ,  Chemistry,  etc),  multi-physics,  multi-scale  approach  to 
complex  space  propulsion  problems 

SUB-AREAS  IN  PORTFOLIO: 

•Coupled  Materials  and  Plasma  Processes  Far  From  Equilibrium 

with  Sayir,  Harrison  (RSA),  and  Luginsland  (RSE) 

•Novel  Energetic  Materials 

Multi-aaencv  Cooidination  Committee:  Petris/DTRA,  Doherty/DHS,  Anthenien/ARO, 
Bedford/ONR,  Spowart,  Hawkins/AFRL,  Palaszewski,  Fletcher,  Sayir/NASA, 
Pagoria/LLNL,  Owrutsky/NRL,  Birkan,  Berman/AFOSR 

•Nonlinear,  multi-scale,  multi-physics  high  pressure  combustion 
dynamics 

with  Fahroo  (RSL),  Darema  (CC),  and  Li  (RSA) 
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Coupled  Materials  and  Plasma  Processes  Far 

From  Equilibrium 

Kick-off  meeting,  NASA  Glenn  RC,  29-30  November  2011 
with  Sayir,  Harrison  (RS A),  and  Lupinsland  (RSE) 


steady-state  powered  low-density  plasmas 

~  1015  #/cm3 

Reverse  field  configuration 


IK  l  i 


Pulsed-powered  relatively  low-density 
plasmas,  ~io13-20  #/cm3 


Plasma/Electrode  Interactions  in  High  Current 
Density  Environments  (HPM  sources)  supercapacitor 

500kV,  10kA,  GW-class  EM  fields,  100ns  to  1  jus 


electric  thrusters 


relativistic 

magnetron 


Pulsed-powered  high-density  thermal  plasmas 

Pressure  ~  1 0’s  M pa,  Temperature  ~1-10eV,  Pulse  time  ~  1  us-1  m 


railgun 


Current  (I) 


/ 


Sliding  Contact  Insulator 


Magnetic  Field  (B):  j  ||  Force  (F)-> 


□ 


IT 


Current  Density  (J)  ■  l*Area 
F  *  J  x  B  =  mass  •  acceleration 


\  / 


Characterize  Surface  /  particle  Interactions  in  Space 
Environment  to  mitigate  contamination,  charging, 
thermal  control,  undesired  optical  backgrounds 


Satellite 

contamination 

/charging 
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steady-state  powered  low-densitv  plasmas 

Understand  and  Model  the  interactions  among  the 
low  pressure  plasma,  material,  and  energy  flow 


TEAM  1 

Nasr  M.  Ghoniem(UCLA) 

Dan  Goebel  (JPL) 

Igor  D.  Kaganovich  (Princeton) 
Yevgeny  Raitses  (Princeton) 
Shahram  Sharafat  (UCLA) 
Brian  Williams  (Ultramet) 
Richard  Wirz  (UCLA) 


What  are  the  relationships  between  surface 
architecture  and  secondary  electron  emission,  and 
the  damage  energy  fluence  limits? 


TEAM  2 

Mitchell  L.  R.  Walker  (Georgia  Institute  of  Technology) 
Alex  Kieckhafer  (Georgia  Institute  of  Technology) 

Jud  Ready  (Georgia  Tech  Research  Institute) 

Greg  Thompson  (University  of  Alabama) 


Howto  model  plasma-material 
interaction  to  characterize  grain 
detachment,  sputtering  leading 
to  plasma  modifications? 


TEAM  3 

Manuel  Martinez-Sanchez  (MIT) 
Mark  Cappelli  (Stanford), 
Dennis  Whyte  (MIT-PSFC) 


What  is  the  the  effect  of  sheath  instabilities,  gas  retention,  and  plasma-induced 
structural  modifications  on  global  performance  ? 
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Multiple  Scales  in  low-density, 
steady-state  powered  plasmas 


c/) 


Sheath 


Thruster 

scale 


Plume  scale 


Timescale 


Electron-neutral 
collision 

Electron-wall 
collision 
Electron 
cyclotron 


Sputtering 

Secondary 

electron 

emission 

A^rrm- 
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Multiple  Scales  in  low-density,  steady-state 
powered  plasmas  -  Computational  Tools 


Length 

scale 


pm  -  mm 
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Multiple  Scales  in  low-density,  steady-state 
powered  plasmas  -  Diagnostic  Tools 
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Timescale 
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Multi-scale  Modeling  of  Materials 
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mechanical 

property 

experiments 
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testing 


micro- 
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B.  Wirth,  2004 
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Example: 

Secondary  Electron  Emission:  Good  or 


•According  to  the  Classic  Fluid  theotx  qL  Hobbs  and  Wesson,  SEE  is  GOOD  //,  leads  to 

reduced  wall  erosion  //  Martinez  -  Sanchez  /  MIT  (1997) 

•Fluid  theory  assumes  isotropic  equilibrium  (Maxwellian), 
and  for  walls  with  very  high  secondary  electron 
emission,  sheath  collapses,  sheat  potential  negligible !! 


Bulk  Plasma 

E= 0 


Wall 


electron  velocity 
distribution  function 

0.2  r 


120 


One  Problem:  can  not 

predict  the  experimental 

results 

• WHY! 


90  ■ 


~  60  - 


E  30 

3 

E 

'x 

& 

0 


Low  SEE  segmented^ 

i  ~ 


o.i  - 


High  SEE  BN  channel 
'fluid  theory  for  high  SEE 


40 


wv  (eV) 


100  200  300  400  500  600  700  800 

Discharge  voltage,  V 

•Plasma  mean  free  path  in  thrusters  are  too  high  to  achieve  equilibrium  (not  enough  electron-electron  collisions) 

•According  to  the  fluid  theory,  maximum  electron  temperature  does  not  change  with  discharge  power  due  to  huge 
electron  heat  flux  to  the  wall 


m 
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w  Example:  Secondary  Electron  Emission:  Good  or  Bad? 

Model  should  account  for  Non-Equilibrium  Effects  !! 


Kinetic  theojx  of  Meezan  and,  Caopelli  /  Stanford.  High  Secondary,  Electmn  Emission 

Depletes  Tail  of  the  Isotropic  Electmn  Velocity.  Distnbution 


•Solution  of  the  Boltzmann  Equation ,  ISOTROPIC 


0  0.1  0.2 


Wall  potential 


Distance  from  exit  (mm) 


electron  velocity 
distribution  function 


w,  (ev) 


20 

0 

wy(eV) 

-20 


•High  Energy  Electrons  lost  at  wall 

•Isotropic,  can  not  predict  secondary  electron  emission  beams,  and  sheath  instabilities 


•Sheath  does  not  collapse,  so  Secondary  Electron  Emission  has  little  effect  on  EROSION! 
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Maximum  electron  temperature,  eV 


t/  Example:  Secondary  Electron  Emission:  Good  or  Bad? 

Particle-ln-Cell  (PIC)  sim  ulations  suggests  that  the 
electron  velocity  distribution  function  is  Anisotropic  !J 

•Correctly  predicts  secondary  electron  emission  beams,  and  sheath  instabilities 
(Raitses,  Kaganovich  /  Princeton) 


120 


90  - 


60  - 


30  - 


Low  SEE  segmented^ 
* " 


o 

electron  velocity 
distribution 
function 


0.1 


0.2 


[inetictfheory 

High  SEE  BN  channel 
fluid  theory  for  high  SEE 


- 1 - 1 - 1 - 1 - 1 - 1 - 

100  200  300  400  500  600  700  800 

Discharge  voltage,  V 


sharp  peaks  are  due 
to  the  Secondary 
Electron  Beams 


40 


Wv  (eV) 


wx  (eV)  '  ,40 

•Secondary  electron  emission  does  not  changethe  sheath  potential 

•Secondary  electron  beams  cause  instabilities  near  the  sheath  surface  due  to  the  “BUNCH  UP” 
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Walls  Made  From  Carbon-Based  Materials  with 
Different  Micro  and  Macro  Structures  Can  Have  Very 
Different  Effect  On  Plasma  and  Sheath  Instabilities 


Diamond  Wall  BHT-200-Cappelli/  Stanford  (Secondary  Electron  Emission  is  unknown) 


•Fundamental  experiments  verified  high  sputter  resistance 
(GOOD!) 

•Diamond  (carbon)  walls  exaggerated  plasmafluctuations- 
olasma  very  UNSTABLE  leading  to  VERYLOWTHRUSTER 

EFFICIENCIES!!! 


Carbon  Velvet  Wall-  Raitses.  Fisch.  Kaganovich  (Princeton)  (has  ZERO  Secondary  Electron  Emission) 


•Carbon  velvet  acts  as  almost  ideal  “black  body”  absorbing  all 
incident  particles  preventing  SEE 

•With  non-emitting  carbon  velvet  walls,  thruster  operates  more 
stable  (no  SEE  induced  instabilities  of  the  plasma-sheath, 
attenuated  breathing  oscillations) 

•With  carbon  velvet  walls,  the  maximum  electric  field  can  be  2-3 
times  larger  than  with  ceramic  walls 


•Same  element,  different  structure  and  architecture 
different  result! 

•Hypothesis:  surface  architecture  affects  performance! 
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•Novel  Energetic  Materials 

Multi-aaencv  Coordination  Committee:  Petris/DTRA,  Doherty/DHS,  Anthenien/ARO, 
Bedford/ONR,  Spowart,  Hawkins/AFRL,  Palaszewski,  Fletcher ;  Sayir/NASA, 
Pagoria/LLNL,  Owrutsky/NRL,  Birkan,  Berman/AFOSR 
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Novel  Energetic  Materials 

Workshop,  23-24  August  2011 ,  Arlington,  VA 


The  Sciences 


Discovery,  understand,  model,  and  exploit 
novel  energetics  materials  to  obtain : 

•  Smart  Responsive  Materials 

•  Nanoenergetics 

•  Energetic  Liquids,  Oxidizers,  and 
Binders 


through  Multiscale  approach  from 

atomistic  to  macroscale  o,deredarrayso, 

nano-porous  silicone 
composites 


& 

Ml 


Next  Step 


graphene  sheets  decorated 
with  energetic  organics  and 
metallic  nanoparticles  for 
performance  enhancement 


encapsulated  nanoscale  fuels 
/catalysts,  nanoporous  fuel  / 
oxidizer  composites  for  control 


surface  functionalization 
particle  morphology,  and 
defect  reduction  to  decrease 
sensitivity 


nPt 

on  graphene 


Current  S&T  Effort 


Frozen  Propellants 
Nickel  Aluminum 
Metal  hydrides 


State-of-the-Art 


Hydrogen  and 
hydrocarbon, 
Ionic  Propellants 
ADN,  HMX 


RDXwith  nano-AI 
encapsulated 


Potential  Impacts 


Polystyrene  coating 
Ammonium  Borane 


•Mission tailored  performance, 
and  burning  rate,  switchable, 
smaller  platforms 

•  Enhanced  and  new  interactions  with 
external  stimuli  for  total  control  of 
reaction 

•  Reduced  sensitivity,  increased  safety, 
and  better  mobility 

•  Enabled  new  missions 


jdfor  public  release;  distribution  is  unlimited 
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High  speed  OH  PLIF  reveals  that  coarse  ammonium 
perchlorate  burns  much  faster  at  high  pressures 


•The  diffusion  flame  structure  changes  from  a  jet-like  to  a  lifted  sheet-like  diffusion  flame  as  pressure  is  increased 
because  of  the  relatively  high  local  burning  rate  of  the  coarse  AP 


1  atm:  Fluorescing  coarse  AP  crystal  is  shown  in  red.  Dashed  line  is  the 


0  ms 


115.6  ms 


165.2  ms 


230  ms 


281  ms 


iSiOSa 

State-of-the-Art  3-D  simulation 
at  high  pressure  (6  atm) 


6  Atm:  The  relatively  fast  burning  crystal  cannot  be  seen  because  it  is  below  the 

surface. 


- 


1  mm 


1  mm 


1  mm 


1  mm 


1  mm 


1  mm 


1  mm 


0.0  ms  1.2  ms  2.0  ms  5.2  ms  65.2  ms  66.0  ms  67.2  ms 

■High  speed  OH  PLIF  also  reveals  that:  Coarse  AP  is  not  affected  by  catalyst  (Fe203  and  CuO)  addition,  Catalyst  should 
be  inside  coarse  AP  in  order  to  have  a  greater  affect  on  performance  4ggp^, 

■  I 
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Nano-scale  features  affecting  meso-scale  behavior 

Example:  mechanically  activated  Aluminum  + 

Fluorocarbon  mixture 


Al  +  polycarbon 
monofluoride 
Powder  Mixture 


Mechanical 

Activation 


Uniform  composition  and  induced 
lattice  defects  (stored  energy) 


S E M/E DX  (Scanning  Electron 
Microscopy/Energy  Dispersive  X-ray 
Spectroscopy) 

Activated  energetic  particles 
Aluminum=blue,  Fluorine=red 
Carbon=green 


25 


DSC  analysis  of 
Al  +  polycarbon  monofluoride 


Al  Fluorination 


20  - 


Lattice  Defect 
Heat  Release 


mechanical  activation 


I 

CD 


Al  Melt^Ai^ 


5  - 


Formation? 

Al  Fluorination  \  \ 


No  mechanical  activation 


-5 


i  i  i  i  i 


100  200  300  400  500  600  700  800 
Temp  (°C) 


•  Fluorinated  graphite  encapsulated  inside  aluminum  at  nanoscales  can  provide 
increased  combustion  efficiency,  reduced  ignition  temperature  and  agglomeration 
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i  ^Ammonia  Borane  (NH3BH3)  as  propellant  additive, 

20%  hydrogen  by  mass,  can  significantly  increase 

rocket  performance 

H 

•  Ammonia  Borane  added  to  hybrid  fuel  (paraffin), 

Up, exp  increased  -10%  with  20%  mass  addition 

♦  Problem:  Significant  AB  addition  led  to  condensed  phase  product  accumulation  on  fuel  grain 

♦MD  simulations,  kinetic  calculations,  and  TGA/DSC  and  Confined  Rapid  ThermoIvsis/FTIR/MS  experiments 


0<SN| 

NH3BH3X 

(BN)x(H)4x 

poly  amino  borane 


hydrogen  eliminations 


+  H- 


oxidation 

V” 

h2o,  h2,  hobo,  n2 


+  02,  OH,  O 


w/o  02  (or  poor  mixing) 
and  /or  lo  w  heating  rates 


Gas  phase  polymerization 
(BN)x(H)4Xi(BN)x(H)2x 


(BN)x(H)2x  l 

polyimino  borane 

Condensed 
phase 
polymerization 

•  dehydrogenation  of  AB  remains  important  in  combustion  processes  at  low  temperature 

•  polymerization  processes  have  lower  energy  barriers  than  dehydrogenation 

•  an  f  in  heating  rate  of  25x’s  results  in  a  jof  2.6x’s  in  mass  accumulation  in  condensed-phase 

♦Hypothesis:  Smaller  ammonia  borane  particles  mav  resolve  the  problem 
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•Nonlinear,  multi-scale,  multi-physics  high  pressure  combustion 
dynamics 

with  Fahroo  (RSL),  Darema  (CC),  and  Li  (RSA) 
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Nonlinear,  multi-scale,  multi-physics  high 
pressure  combustion  dynamics 

Workshop,  15  June  2011,  London;  Workshop,  23  August  2011,  Arlington,  VA 


PARADIGM  SHIFT  IN  SIMULATION: 

•Reduced  Basis  and  Stochastic  Modeling  of  a  Complex  High 
Pressure  Combustion  System  to  identify  physical  mechanisms 
responsible  for  the  observed  dynamical  behavior 


PARADIGM  SHIFT  IN  VALIDATION: 


Closed-loop  actively  controlled  real-time  hybrid  approach 


Experimental  domain 
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waves  structure  in  a  high  pressure  toroidal 
cavity  with  concentrated  heat  release  zones 

Zinn,  Yang,  Neumeier /Georgia  Tech,  and  Law  (Princeton) 


•Consider  a  toroidal  cavity  with  a  single  combustion  point  source.  The  acoustics  is  governed  by  the 
homogeneous  wave  equation  except  in  the  singular  point. 


•Consider  now  the  homogenous  wave  equation: 

dPa/dt2-C28P'2/dx2  =  0 


XR  -x-ct  Right  propagating  characteristics 


XL  -x+ct  Left  propagating  characteristics 


Right  going  pressure  wave  pR=f(AR) 


Left  going  pressure  wave  pL  =  gUL) 


/(/) 


u  = 


u 


pc 

.L  _  gqL) 

pc 


•With  combustion: 


The  compact  heat  release  source 
modifies  the  out  going  waves 


/l=/'+zF?-*' 


SW  2o 
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^  k  The  increased  mean  flow  due  to  the  convectinq 

^^40  standing  waves  causes  spinning  waves  with  increasing 

amplitude  leading  to  instability 


Initial  conditions  If  no  tangential  mean  flow:  spinning  to  a  Standing  wave 
Right  travelling  wave  f 


Initial  conditions 
Left  travelling  wave  g 


-04  -02 


0.2  04 


TANGENTIAL  MEAN  FLOW  due  to  recirculation  near  the  injectors  and  inhomogeneous  distribution  of  heat  release,  leading 


Initial  conditions 
Right  travelling  wave  f 


to  convectinq  to  a  spinning  wave 


Initial  conditions 
Left  travelling  wave  g 


•Explained  experimental  results  obtained  at  NASA  (Marcus  Heidmann  ,  1969) 

Zinn,  Yang,  Neumeier  /Georgia  Tech  DISTRIBUTION  A:  Approved  for  public  release;  distribution  is  unlimited 
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Summary 


Space  Propulsion  Portfolio  has  become  a  platform  for  Multi-disciplinary 

activities  in  all  scales: 

•Propulsion 
•Materials 
•Interface  Sciences 

•Plasma  and  Electro-Energetic  Physics 
•Applied  and  Computational  Mathematics 
•Chemistry 
•and  others... 


Will  provide  the  scientific  foundation  to: 

•Reduce  Fuel  Demand  in  Space  /  more  efficient  power  generation  and 
energy  utilization,  increase  spacecraft  lifetime,  and  reduce  /  control  waste 
heat  /  provide  novel  design  propellants  /  increase  reliability  and 
performance 

•understand  and  manage  High  Energy  Density 


•Enable  high-energy  storage  in  ultra-capacitors  with  nanostructured 
components, 


•and  others... 
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